arXiv:nucl-th/9701014vl 7 Jan 1997 


CONFINEMENT AND SCALING IN DEEP INELASTIC 


PROCESSES * 

S.A. Gurvitz 

Department of Particle Physics, Weizmann Institute of Science, Rehovot 76100, Israel 


The existing data for hadron structure functions, F 2 {x,Q^), show considerable Q^- 
dependence, which is mainly attributed to the QCD logarithmic corrections to Bjorken 
scaling. However, at a; ^ 1 the scaling violations are dominated by power corrections 
oc 1/Q^ (higher twist and target mass effects): 

F2(x,Q^-)^Fr(x,Q-^) + ll^ + ---, ( 1 ) 

where F 2 ‘^{x,Q‘^) = F 2 {x,Q‘^ 3> |H(x)|) and the remaining Q^-dependence in F 2 ^(x,Q^) is 
to be attributed to QCD logarithmic corrections only. 

Power corrections can be incorporated in the first term of Eq. (|l|) by using a different 
scaling variable, 

x = (j){x,Q‘^)=x+^^ + ---, ( 2 ) 

so that 

F^{x,Q^) = Fj (r ■(*,(?"),(?") =: Fr(i,Q"). (3) 

The coefficient B, which determines the value of the power correction in Eq. (|I]), is thus 
related to the structure function by B{x) = b{x) Q'^)/dx. In fact, an analysis of data 

in terms of an appropriate scaling variable appears to be more convenient, than the direct 
evaluation of power corrections. 


*WIS-97/2/Jan.-PH. Invited paper given at ISHEPP, Dubna, September 1996. 


1 





In general, the power corrections are generated by confining interaction of partons in the 
final state. At first sight, such an interaction should influence the structure function very 
drastically. Consider for instance the example of two nonrelativistic “quarks” of mass m 
interacting via a harmonic oscillator potential []^ . These quarks are never free and therefore 
the system in the final state possesses a discrete spectrum. As a result the structure function, 
F{q, I/), as a function of the energy transfer z/, is given by a sum of 5 - functions. Obviously, 
it looks very different from the structure function obtained in the impulse approximation, 
which considers the struck parton as a free particle in the final state. This paradox can be 
resolved by introducing a (nonrelativistic) scaling variable y 



( 4 ) 


Then expanding the structure function F{q,y) = F{q, z/) in powers of 1/g, one finds in the 
limit g —cx) and y =const that the 5-peaks merge to a smooth curve, F{q,y) ^o{y), 
which coincides with a free parton response |^. Although this result appears to confirm 
the parton model picture, it does not imply that the interaction in the final state is not 
important. The latter has been merely incorporated in Fo{y) by an appropriate choice of 
the scaling variable y, as shown by Eqs. m- The remaining contribution from higher-order 
(~ 1/q) terms are thus minimized. However, a non-optimal choice of the scaling variable 
could result in very large or even singular corrections to the structure function. One can 
anticipate that an appropriate choice of of the scaling variable is especially relevant at large 
X, where lower-lying excitations should play an important role. 

A general analysis performed in the framework of Bethe-Salpeter equation shows that in 
an analogy with the nonrelativistic case the higher twist terms from a local confining 
final state interaction and target mass effects can be effectively accounted in the modified 
the struck quark propagator, where the quark mass has the same off-shell value before and 
after the virtual photon absorption [Q, 
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FIGURES 



FIG. 1. Diagrammatic representation of the leading contributions to the structure function 
using light-cone variables. Quarks and gluons are shown by solid and wavy lines respectively. The 
modified propagators are marked by 

As a result, the Bjorken scaling variable x is replaced by a new scaling variable x = 
x{x,Q‘^), which is the light-cone fraction of the off-shell struck quark. Explicitly, 

X -\- — y(1 — xY IQ^ 

1 + + 4M2a;Vg2 

where M is the target mass and rus is the invariant mass of spectator partons (quarks and 
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gluons). For ^ oo or for x —0 the variable x coincides with the Bjorken variable x. 
However, at hnite these variables are quite different. 

It follows from Eq. (^) that x depends on the invariant spectator mass, rris- In terms of 
light-cone variables. Fig. 1, it can be written as 

ip'± -P±f 


mj = mQ + {1 — z)- 


z' — z 


+ - z) . , + P i - pi. 


( 6 ) 


where mo is the diquark mass, and 1 > z' > z. We approximate m^ as an effective spectator 
mass depends only on external momenta. Since ^ x ~ x and z' ~ one gets from Eq. 


m^ ~ mgC(x, (5^)(1 — x) (7) 

Eq. d^) for the invariant spectator mass looks quite appealing apart from its relation 
to Eq. (I). Indeed, x = 1 corresponds to elastic scattering, when no gluons are emitted. 
Therefore in this case the spectator is represented by a diquark. When x decreases, gluons 
are emitted and m^ increases oc (1 — x). The coefficient C'(x, in Eq. determines the 
rate of increase of the spectator mass with and x. It can be found self-consistently from 
the evolution equation. However, when x ~ 1, one can take C(x,(5^) — C'(1,Q^) ~ const, 
because of Q^-dependence of the spectator mass is less important than its x-dependence 
near the elastic threshold. We roughly estimated the value of C by using the Weizsacker- 
Williams or “equivalent photon” approximation, utilized in Ref. for derivation of the 
evolution equation. One hnds from ^ that the light-cone fraction of the “equivalent” gluon, 
z — z', (Fig. 1) is of order asln((5^/Qo) region of large x. However, the probability 

of the gluon emission is also about the same order of magnitude. Then, as follows from Eq. 
(^), C ~ {{p'±_ — so that one could expect to hnd C on the scale of (GeV)^. In the 

following we regard it as a phenomenological parameter, determined from the data. 

Let us consider the nucleon structure functions in the region of large x, where the power 
corrections to the scaling are dominant. At present, the only available large-x data for 
proton and deuteron structure functions, F^x, are the SLAG data 
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taken at moderate values of momentum transfer, < 30 (GeV/c)^. (The nucleon structure 
functions for higher values of momentum transfer < 230 (GeV/c)^) are extracted from 
BGMDS 0 and NMG data, yet only for x < 0.75). The SLAG data for the proton and 
deuteron structure functions for x > 0.7 and 5 < < 30 (GeV/c)^ are shown in Fig. 2 as a 

function of x. Also shown is the value of Tf (x, and F^ix, for Q^= 230 (GeV/c)^ and 
X = 0.75 taken from the BGDMS data |P]. The data points close to the region of resonances 
were excluded by a requirement that the invariant mass of the final state (M + 
is greater than (M + A)^, where A = 300 MeV. In addition, we excluded the data points 
with X > 0.9 from the deuteron structure only (Fig. 2b). The reason is that the deuteron 
structure function can not be represented as an average of the proton and neutron structure 


functions for x > 0.9. Indeed, the calculations of Melnitchouk et al. show that the ratio 
2F2/{F2 + F 2 ) is about 1.13 for x = 0.9 and Q^=5 (GeV/c)^, and it rapidly increases for 
X > 0.9. However, for x < 0.85, this ratio is within 5% of unity ||I2[ . 


One hnds from Fig.2 that the structure functions show no scaling in the Bjorken variable 
X. Also, very poor scaling is obtained when the data are plotted as a function of the 
Nachtmann variable ^ [^. However, the situation is different if we display the same data as 
a function of the variable x, Eq. (^. It appears that the scaling in the x-variable is strongly 
dependent on the value of diquark mass, mo, Eq. 0 . but is much less sensitive to variation 
of the coefficient C. For instance, the data as a function of x display very poor scaling for 
mo > 600 MeV, i.e. by considering the spectator as build up from constituent quarks. 
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FIG. 2. The SLAG data[5,6,7,8] (5 < < 30 (GeV/c)^) for proton (a) and deuteron (b), are 

shown as a function of the Bjorken variable x. Three high-statistics data sets[7] for ~5.7, 7.6, 
and 9.5 (GeV/c)^ are marked by “x”, and respectively. The point at = 230 (GeV/c)^ 
and X = 0.75 is from ref. [9]. 
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FIG. 3. The data of Fig. 2 are shown as function of the x(x, Q ‘^)—the scaling variable of Eq. (5) 
— assuming mo = 0 MeV and C = 3(GeV)^ for the spectator mass m^, Eq. (7). 

On the other hand, the scaling is very good both for the proton and deuteron data, by 
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taking uiq = 0, i.e. by considering the spectator build up by current quarks [Q]. The results 
are shown in Fig. 3, where the data are plotted as a function of x for mo = 0 and C = 
3 (GeV)^. Note, that the high-Q^ data points from BCDMS data P are very close to the 
SLAG data points, taken at much lower values of Q^. 

Now by using the scaling variable x, Eq. (^) for tuq = 0 we can analyze the nucleon 
structure functions for smaller values of x, where both power and logarithmic corrections to 
the Bjorken scaling play an important role. In this region (0.35 < x < 0.75) the existing 
BGDMS i and NMG m data are extended up to much larger values of momentum transfer 
than the previously considered high-x SLAG data. It allows us to check our predictions in a 
wide range. The results are shown in Figs. 4 and 5 for proton and deuteron structure 
functions respectively. The data points are from SLAG and BGDMS data bins The 
dotted lines show the Q^-dependence of the structure functions due to power corrections 
only. 






FIG. 4. The proton structure function F 2 {x,Q‘^) is shown as a function of at different 
X-values. The dotted lines include power corrections only. They are evaluated according to and 
the scaling variable x of Eqs. (5), (7) with mo = 0 and C = 3 (GeV)^. The additional QCD 
logarithmic corrections evaluated at NLO for different A scales are shown by the dashed (A = 100 
MeV) and continuous lines (A = 200 MeV). 

The total Q^-dependence of structure functions due to the power and the logarithmic 
NLO corrections, is shown by the dashed and continuous lines for A = 100 MeV and A = 
200 MeV respectively. The QCD (logarithmic) evolution corrections are taken into account 
at Next-to-Leading Order (NLO) evolving back, in Q‘^, the structure functions starting from 
an asymptotic value of momentum transfer where the condition F 2 {x, Q'^) ~ Q‘^) (cf. 

Eq. (|lD) is fulhlled (in the present case we choose = 230 (GeV/c)^, which is the highest 
value of the momentum transfer in the BCDMS data |P). 
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FIG. 5. As in Fig. 4 for the deuteron structure function 


One finds from Figs. 4 and 5 that the experimental data are reproduced in a large Q^- 
range for both values of A, although the agreement is slightly better for A = 100 MeV. In 
addition, since the results are strongly dependent on the spectator mass (|^, it is remarkable 
that the same parameters mg = 0 and 0 = 3 (GeV)^ do reproduce the Q^-behavior of the 
structure functions both for large and moderate x-values. 

The above analysis of the nucleon structure functions allows us to extract from data 
the previously unattainable information on the asymtotic ratio of F^(x)/Ff’(x) at large 
X, directly related to u and d-quarks distribution. Indeed, according to the quark-parton 
model, the structure functions of hadrons in the Bjorken limit > cx) and 

X = Q^/2Mz/=const) are directly related to the parton distributions qi(x). For instance 


F2(x,Q^) ^ F2(x) = ^e-xqi(x). 


( 8 ) 


where the sum is over the partons, whose charges are e*. In the region x —*• 1, the contribution 
of sea quarks can be neglected. Then assuming the same distribution of the valence quarks, 
one easily hnds from Eq. that the neutron-to-proton ratio F^^x) /F^i^x) approaches 2/3 
for X —1. If, however, the quark distributions are different, one can establish only upper 
and lower limits for this ratio, 1/4 < F^/Ff < 4, which follow from isospin invariance 


In order to check the parton model predictions one needs to take the structure functions at 
high since the higher-twist corrections to the scaling are very important at high-x region. 
At present, high structure functions — 250 (GeV/c)^) extracted from BGDMS 0 
and NMG data are available only for x < 0.7. The ratio F^ix)/F 2 {x) obtained from 


an analysis of these data shows steady decrease with x. Thus, it is usually assumed 
that this ratio would reach its lower bound, F^/F|’ —> 1/4, for x ^ 1 (although the recent 
analysis []^ suggests that this ratio would be larger than 1/4 for x —> 1). This corresponds 
to d{x)/u{x) —0 for X —1, where d{x) and u{x) are the distribution functions for up and 
down quarks in the proton. 
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Since our scaling variable x takes effectively into account the higher-twist corrections, it 
allows us to reach the asymptotic limit already at moderate values of momentum transfer, 
< 30 (GeV/c)^. Thus, Ff and F 2 extracted from the SLAG data at large values of x P-§|, 
Fig. 3, would represent the asymptotic structure functions: F 2 {x) = F 2 {x) = F 2 {x,Q^ —>• 
cx)), which can be used as an input to determine the ratio F^^x)/F^^x). For this purpose we 
parametrize the asymptotic structure functions as Ff= exp(— J2i=o cbiX^) and determine 
the parameters ai from the best £t to the data in Figs. 3a, 3b. The resulting F^^/Ff ratio is 
shown in Fig. 6 by the solid line. The dotted lines are the error bars on the fit, which combine 
statistical and systematic uncertainties. The dashed line corresponds to F^/Ff = 2/3. For 
a comparison, we show by the dot-dashed line a polynomial extrapolation of this ratio to 
large x, obtained from BGDMS and NMG data by assuming that F^/Ff —1/4 for x —> 1 
H- structure functions 
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FIG. 6. Neutron-to-proton structure function ratio at large x. 


Our results shown in Fig. 6 demonstrate that contrary to earlier expectations, the ratio 
does not approach its lower bound, but increases up to approximately 2/3. The 
latter is the quark model prediction, assuming the same distributions for each of the valence 
quarks. The accuracy of our results will be checked in future experiments, which will provide 
high data for the structure functions at large x. 

The described above method for evaluation of the higher-twist corrections using modihed 
scaling variable can be very useful for not only for an analysis of high-x nucleon structure 
functions, but also for a treatment of many different problems, where the interaction of 
partons in the hnal state becomes important. In particular, I would like to mention the 
sum rules, where the role of higher-twist effects remains an open problem, and the structure 
functions at x —0. In this region the spectator mass is not described by Eq. (0), valid 
for large x. If oc 1/Q^ for x —>• 0, then (x — x)/x oc 1/Q^, Eq. and the role of 
higher-twist correction would be substantial in this region too. 
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